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SYNOPSIS 

The thesis entitled "SIMULATION OF SURFACE 
GENERATED IN FINE GRINDING - A PROBABILISTIC APPROACH" 
deals with the generation of theoretical surface profile 
simulating the actual workpiece surface profile obtained 
under single pass plunge-cut grinding conditions through 
estimation of probabilistic parameters to describe the 
cutting, .surf ace. A theoretical model is proposed to 
describe the probability density functions of the random 
positions of grains on the wheel surface based on the 
interference of grooves produced by different grains, 
effective profile depth, and the elastic grain deflection 
The workpiece surface being the sample realisation of 
random grinding process, power spectral density of the 
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simulated and experimental profiles are compared which 
provides information about the real system in terms of 
radial distribution parameter and effective profile depth 
for various cutting conditions. The surface roughness 
and active grain density are also calculated as outputs 
of the grinding system. 

Chapter I introduces the earlier developments 
of both the deterministic and probabilistic analyses of 
the grinding process and outlines the need for the pre- 
sent investigation. 

The model based on the probability density func- 
tion of the random position of grains on the wheel sur- 
face, and the generation of the theoretical profiles in 
the longitudinal as well as in the transverse direction, 
for various grinding conditions, based on the interference 
of cuts made by random grain tips is described in 
Chapter II. 

Chapter III provides the statistical analysis 
to obtain the functional dependence of radial distribu- 
tion parameter and effective profile depth on the table 
speed and -wheel depth of cut, by comparing the power 
spectra of the simulated and experimental profiles using 
a systematic variation of the table speed and depth of 
cut through a 2 x 2 factorial experiment. 
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Numerical techniques to evaluate the power 
spectral density and generation of surface profile with 
all random possibilities of interaction of grooves is 
discussed in Chapter IV. Digitization and handling of 
random data uging standard methods are also described 
inlthis chapter. 

Chapter V presents an approximate analysis to 
obtain the surface roughness and variation of active 
grain density with undeformed chip-thickness. 

Experimental details and results are summarised 
in Chapter VI. It is shown that the effective profile 
depth is independent of depth of cut and table speed, and 
that the radial distribution parameter is linearly rela- 
ted to table speed and depth of cut. Results indicate 
that the power spectra of the simulated profiles obtained 
using the proposed model are in close agreement with the 
experimental profiles, and shows the dependence of sua>- 
face roughness on the cutting conditions. 



CHAPTER I 


INTRODUCTION 


1 .1 Introduction 

In recent years, grinding has received more 
attention because of the ever increasing higher trend 
towards precision in the processing of varied and stronger 
materials which are difficult to machine. Increased inte- 
rest in the mechanics of grinding as a metal cutting 
process has contributed significantly in understanding its 
fundamental parameters. Though, the grinding process, in 
many ways, is similar to milling, it is unique since the 
material removal is carried out by small closely spaced 
and randomly placed grains of hard abrasive materials. 

The art of using abrasives is very old and is 
being successfully practiced since the pre-historic times. 
Archaeological evidences confirm that early man in the 
paleolithic age, from 25,000 B.C. to 15,000 B.C., was well 
versed in the use of abrasives in fashioning and sharpen- 
ing of bones and horns pi J . Many synthetic grains such 
as Aluminium Oxide (produced by fusing Bauxite in a 
special electric furnace |2{), White Aluminium Oxide (a 
more refined form of Aluminium Oxide), or Silicon Carbide 
have Replaced the use of natural abrasives such as emery, 
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and corrundum. Aluminium Oxide is recommended for grinding 
metals of high tensile strength due to its toughness. 

Hhite Aluminium Oxide has the characteristics of high 
friability and first rate cool-cutting properties 3 } • 

Silicon Carbide , the first synthetic abrasive discovered 
by Archeson around 1395 [5 j ? was at one time the hardest 
manufactured abrasive but it is very brittle j~6, 7, 8_ . 

More research is continuing into the development of new 
types of abrasives. 

Abrasive grains are held together by bonding 
materials. Majority of the wheels are made with vitrified 
bonds. Resenoid bonded wheels are extensively used on 
special swing frame, floor stand, and high speed portable 
grinders for the grinding of steel, malleable and iron 
castings where quick removal of metal is the primary 
consideration. 

Grinding wheels are made in a variety of grain 
size, grade and structure. The important properties of 
an abrasive grain (or grit) are hardness | body strength, 
or toughness (inversely proportional to friability), and 
attrition resistance jjO, 1 TJ . Hardness values are ex- 
pressed by Woodddll scale pi 2 ] which are directly compara- 
ble with cohesive energy" density pi 3j * Grade is the 
degree of hardness of the wheel. It varies according to 
the tenacity with which the bond holds the abrasive grits 
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together and does not represent the hardness of the 
abrasive material itself. The grain size is designated 
in terms of openings per linear length in the screen to 
size the grain. The sizes usually range from 10 to 120; 
10 being very coarse and 120 being very fine. The struc- 
ture of the grinding wheel denotes the spacing of grain 
and controls the density of the wheel. The grinding 
wheel is not a dense mass of abrasives and bond. It has 
intergranular space distributed throughout its structure 
which helps to clear the wheel face of metal and remove 
the chips cut by abrasive grains. 

The grinding process differs fundamentally from 
other cutting methods in the following j~lk 

(i) the space distribution of the cutting edges on 
the cutting surface of a grinding wheel is of 
random nature; 

(ii) the shapes, sizes, and orientation of cutting 
edges vary over a wide range; 

(iii) the radius of curvature of cutting edge of a 
grain determines chip thickness; 

(iv) grains, held elastically by bond bridges are 
displaced in tangential and normal directions 
during grinding; 




if 

(v) the types of wear of the cutting edges are 
quite different from those of other cutting 
tools, 

(vi) the chip thickness is very small with wide 
variations in shape and size; 

(vii) the temperature of cutting edges are extremely 
high. 

Because of the large number of parameters involved, 

an enormous number of trials will be required to determine 

*» 

experimentally the relative importance of each parameter 
influencing the grinding process. Efforts have been made 
to obtain a relevant basic parameter to combine several 
effects produced by a large number of parameters. The 
various proposals related to basic grinding parameters are 
split into two groups; one based on geometric relations 
(speed, depth of cut, cutting edge spacings, etc), and the 
other based on concepts such as specific energy, specific 
contact pressure etc. \v? - 23 • It is generally agreed 

that the basic parameter should be decided in such a manner 
that it refers explicitly to those values which can be 
evaluated with reasonable accuracy. The parameter unde- 
formed chip thickness, t, is generally accepted as the 
basic parameter. 
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Some typical kinematic arrangements used in 
grinding such as (a) internal grinding, (t>) horizontal 
spindle surface grinding, (c) external cylindrical grind- 
ing, and (d) abrasive • cut-off are shown schematically in 
Figure 1.1. 

1.2 The Role of Chip-thickne ss 

The important factor in deriving the expression 
for undeformed chip-thickness, t, is the number of cutting 
elements per unit area of wheel surface, C, that are pre- 
sent in a layer whose thickness is related to that of 
thickest chip produced [ 2*+-27 j . The deeper the grain 
lies above the outermost wheel surface, the greater is 
the probability that a grain will not contact the work- 
piece surface because a more protruding grain immediately 
ahead of it has already removed the material. The other 

parameter, affecting t is the ratio of mean scratch width 

*\ 

to scratch depth, r, { 28 ~j. The chip size determines the 
area of contact between chip and the grain and the force 
on the entire grain. The specific energy, U 
given by 

F» ¥ 

0 = ( 1 . 1 ) 

vd 

f 

where F is tangential force per unit grinding width, 

3 ? ... 

V and v are the wheel and table sueeds re selectively and 
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d is wheel depth of cut. Figure 1.2 shows a plot of 
U versus t on a log-log scale £2*fj. In plotting this 
curve the values of t were computed using equation 
(1.10) in which C was assumed to have a constant value 
depending on the grain size. 

■ t ■ ■ ■ 

1*3 Mechanics of Grinding 

As the average linear grain spacing, f\ , between 
two grains on the wheel surface in the same peripheral 
plane passing through their tips can not be measured 
directly, equation for t can be obtained from volume 
continuity £ 28 ] . In surface grinding, represented schema- 
tically in Figure 1.3 a, the actual shape of the chip 
produced is shown in Figure 1 .3 b* An idealised chip 
shape (Figure 1.3 b) can be assumed to ensure volume con- 
tinuity. Since the chip length, l c , is much greater than 
t and s f where is feed per grain, given by 

Sf = AvA .(1.2) 

vrj , the volume of the chip shown in Figure 1 .3 c ifis given 
approximately by 

v.| = 2 l c ^ k 0*3) 

where b is the average width of the chip. The volume 
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removal rate by the number of cuts made per unit time. 


_ jSLJlja _ 2_d 
“ V b C “VC 


<1.40 


where C is the number of active grains per unit area 
on the wheel surface. From equations (1.3) and (1.4) 
we get, 


2 v d 


V C 1 b 
c 


0 . 5 ) 


Ignoring the local wheel deflections, 


l c - VB e d 


( 1 . 6 ) 


where D is the eauivalent wheel diameter given by 
6 


D D 


D 


w 


D + D 


(1.7) 


w 


In equation (1.7), D is the wheel diameter, D w is the work- 
piece diameter, positive sign is for external and negative 
sign is for internal grinding. For plunge cut surface 


grinding, D is infinite. Therefore, 


D d 
e 


Hence, 




" V C bVD 
Assuming [ 2 2, 24, 29 


(1 * 8 ) 


O .9) 
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t 

f 



( 1 . 10 ) 


Expressing C in terms of .average linear grain spacing. 


t, can be obtained from equation ( 1 . 9 ) as 

t = 2 A CvA)\/|- (1.12) 

In the above analysis, the path traced by an individual 
grain has been assumed to be circular instead of trochoi- 
dal. This assumption seems to be justified in view of 
the speed ratio, q = V/v, being very high. In equation 
'(1.10), the unknown parameters are C and r. 


1 Determination of C a nd r 

The number of cutting edges per unit length or 
unit area, the distance between cutting edges, and the 
distribution of cutting edges in the direction normal to 
the wheel surface quantitatively characterise the cutting 
surface, C was originally determined by rolling a dressed 
wheel over a soot-coated plate and counting the number of 
contacts indicated on an enlargement using the plate *as 
negative^ 2*f, 25~j . This method gives a very high value of 
C as vibrations, local wheel and work deflections, etc., 
are not included. Several other methods have been used 

0, 3l] to observe 


to evaluate C, It has been suggested 13 
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the wear flats produced on the grains and to get the C 
value by dividing the number of wear flats by the viewing 
area, or by counting the number of chips produced after 
separating the chips from the debris collected J^3 2j * The 
latter method is based on the assumption that each cutting 
element produces only one chip. In another method |33| , 

i— ~ -I 

C was evaluated from the length and spacing of the scrat- 
ches produced during grinding on a lapped inclined surface. 
A plot of C versus z‘ (Figure 1.4) has been obtained, where 
z' is the distance measured radially inward from the outer- 
most grain on the wheel face. From Figure 1.4, z 1 is taken 
equal to t, and the value of C is obtained. The equation 
(1.10) is solved iteratively till the assumed and calcu- 
lated values of t match. 


It has been observed j 34 j that when the value 
of U is evaluated for any value of t obtained in con- 
junction with Figure 1 .4, the plot of U versus t 
(Figure 1.5) has a very steep slope, which suggests that 
there is practically no difference in t for the finest 
and ccfefibsest grinding . The difficulty is believed to lie 
in the fact that the curves of Figure 1.4 have been extra- 
polated in obtaining Figure 1.5. Therefore, it has been 
suggested [~ 29 ~] that the curves of Figure 1.4 will not 
continue upward indefinitely as shown by dotted curve in 
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Figure 1 .6 but they should level off beyond some value 
of z’ (or t) as shown by dashed curve. Thus, in the 
practical range of values of t, C will be approximately 
constant. Attempts have also been made to find the re- 
lationship between the size distribution of abrasive 
particles removed from grinding wheel by dressing and the 
number of active grains per unit area of the wheel for 
different grain sizes and wheel grades j 3 

The quantity, r, was originally obtained from 
a taper section j_2Vj . It depends on the shape of a typi- 
cal grit which has been the subject of much uncertainty. 
Earlier, it had been a widely held belief that abrasive 
grits generally have a large negative rake angle [ 36j. 
However, this concept cannot explain a number of impor- 
tant features of the grinding process such as the exis- 
tence of large specific energy for metal removal. The 
action of a single abrasive grain were simulated by many 
researchers using spherical and conical tools L37“ 1+ 3J* 
However, the manner in which tools having these idealised 
shapes form a groove or remove a chip is still not clear. 
These workers jjJ7-^3j| suggest that grinding mechanism can 
be broadly divided in two categories, fine and rough 
grinding. Hence, if grinding is viewed as a scraping 
process, two typical grit geometries can be assumed. At 
large wheel depths, the grit may be assumed to be conical 
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•whereas for small wheel depths, a spherical shape may he 
a better approximation j~44j . Tests conducted with single 
grains j 45, 46~j indicate that in fine grinding, the average 
grit shape may be approximated by a sphere. Based on this, 
the grinding mechanism is assumed to be of extrusion type 
rather than concentrated shearing mechanism involved in 
normal metal cutting operation P+zjJ • It has been shown 
that scratches produced in overcut fly-milling is closely 


approximated by an arc of a circle [_48j , and the radius 

of curvature is independent of grain type and grinding 

conditions in fine grinding but varies with the grain size 

(Figure 1.7)* The equation for undeformed chip thickness 

- - 

can be rederived as 1 48 


t = ■■■ ( 1 . 13 ) 

V C \/D V2> - d 

■where p is the grain tip radius as shown in Figure 1.8. 
Single grain abrasive studies suggest that the best esti- 
mate for the average transverse grain shape is a cone with 
radius jP at the tip, thus for material removal, the 
conical shape is assumed but for surface finish analysis, 
the spherical shape is assumed f 49 j. 

1 .5 Grinding as a Random Process 

In most of the material removal processes, a 
single point or multipoint cutting tool with well defined 
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cutting edge geometry is used. For grinding process, 
unfortunately the cutting edges are not well defined or 
equispaced, rather they are placed randomly on the grind- 
ing wheel. Thus they produce a surface whose behaviour 
can not be predicted exactly in a deterministic manner. 

Due to this randomness of cutting particles on the grind- 
ing wheel surface, it is natural to strive for a stochas- 
tic model of the grinding wheel and the grinding process. 

1 .6 Stochastic Model 

In most of the engineering problems involving 
randomness, the first step is to consider randomness in 
the analysis and ultimately to construct a reasonable 
stochastic model representing such randomness which gives 
the output characteristics closest to desired output. It 
has to be reasonable in the sense that it can be handled 
analytically without too much cost of computation and should 
also be consistent with the quality and quantity of infor- 
mation pertinent to the particular random phenomenon of 
interest. An overly sophisticated and complex model that 
cannot be verified by means of existing data and possibly 
by way of sound engineering judgement could simply compli- 
cate rather than solve the problem at hand [~5oJ« Stochasr- 
tic modelling of physical system can be achieved through 
Monte-Carlo Simulation method or Markov-Chain method. 
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Monte-Carlo method j 51 j is based on random 
numbers generated for a known or assumed probability 
density function of any variable. Depending upon the 
dynamics of the process, the random response of the sys- 
tem is evaluated and compared with desired output of the 
system. 

Markov process Q>2~] is best visualized as a 
process where the stochastic dependence of the future is 
independent of the past information given the present 
state precisely. If the state of space is discrete or 
contains finite or countably infinite number of states, 
the process is known as Markov- Chain. 

1 .7 Previous Models 

Surface roughness of the workpiece is one of 
the most important output parameters in precision machi- 
ning, The grinding process can therefore be analysed in 
terms of the surface roughness produced on workpieces. 
Yang and Shaw have derived an expression for centre 

line average (CLA) value of the surface roughness for 
surface grinding, assuming that all the grains are at the 
same radial position on the grinding wheel generating 
chips of equal size. The mean peak to valley height has 
been evaluated by assuming either the grain population 
density to vaiy parabolically with the distance from the 
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outermost grain on the wheel face j" 54 j , or the grains 
to he uniformly spaced, j 33 J • The surface roughness values, 
thus, obtained 53 s 54 , 33 J were either too low or high ...as 
compared with the actual measurements. This was primarily 
because these analyses were based on the assumption of 
some form of population density function of grain position 
with radial distance without considering the dynamic effect 
of table speed, v, wheel depth of cut, d, and elastic 
grain deflection. The grains responsible for surface 
generation are distributed within an effective profile 
depth, d , which is defined as the radial zone on the wheel 
surface within which the grains contribute cutting 55j« 
Comparison of surface roughness in terms of CLA value or 
mean peak to valley height is not likely to give a good 
representation of the surface generated, since the sur- 
face profile itself is not considered. In this respect, 
autocorrelation and Fourier spectrum analyses are likely 
to lead to a much better representation of the generated 
surface profile. This can be obtained through a random 
model . 

1 .8 Random Models 

A simple random model was proposed by Baul 0>£j 
where the asperities and .valleys of abrasive surface have 
been simulated by rods of varying heights having a 
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truncated normal distribution. The sample distributions 
of the distances between active asperities, defined as 
grit pitch, and the depth of penetration were then obtain- 
ed. The cutting performance in belt grinding was studied 
theoretically 57 j t>y analysing the profile using Markov- 
Chain theory. The theoretical distribution for distan- 
ces between two active grains was obtained by observing 
the lands and voids of the representative profile of the 
abrasive surface through the transition probability matrix 
p of going from state o to 1 . When the profile height 

o 1 

measured at constant intervals along some fixed length 
exceeded a fixed height, it was defined as state 1 and 
when it failed to exceed, was state o. 


The grinding process was simulated by Yoshikawa 
and Sata j 58~1 using the Monte-Carlo method. The transfer 


function of the process was obtained using the simulated 
profile as the input and the experimental profile as the 
output. Figure 1.9 shows that the transfer function of 
the process, thus obtained, is very much amplified at 
lower frequencies and almost unity at higher frequencies. 
The reasons attributed for the distortion of transfer 
function at low frequency and the assumption that the 
radial distribution of grain tips is proportional to the 
depth from the outermost grain on the wheel face needs 
justification. 
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Peklenik | 59 investigated, the grinding process 
as a linear system -with random input and output. The 
elementary cutting profile of the grinding wheel was taken 
to he a stationary stochastic process described in terms 
of averages and auto-correlation functions. These charac-. 
teristics were obtained for representative grinding wheels 
and were functions of physical and geometrical behaviour, 
and the dressing conditions. The frequency transfer func- 
tion of the surface roughness and the wear of grinding 
wheel provide an insight into the tranf ormation of the 
disperson spectra in the grinding process. Information 
on the shape of transfer function make it possible to 
determine the cutting action of the elementary cutting 
profile in terms of the distribution of the instantaneous 
chip size. 

Autoregressive moving average model was developed 
by Deutsch and Wu j 60~] to represent grinding wheel profiles 
for different combinations of sampling parameters including 
the sample interval, the number of observations, and the 
length of record. The dioice of sample interval determines 
the appropriate model form, the model order and the esti- 
mated parameters. The use of larger number of observations 
is necessary only to increase the confidence associated 
with the estimated parameters. 
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Another simulated model was proposed by Law 
and Mu [ 61~] by systematically varying the table speed, 
wheel depth of cut, and grain apex angle. In this analy- 
sis, the dynamic influence of table speed, depth.' of cut 
over a preselected region and elastic grain deflections 
were predicted through the radial distribution parameter. 
It was proposed that the radial distribution parameter 
has two range of values, one less than 8 for fine dress- 
ing and the other greater than 8 for coarse dressing. 

\ 

This was explained in terms of change in the grain apex 
angle due to dressing conditions. 

In all the previous analyses 55-60 i , the 
transverse shape of the grain was assumed to be triangular 
with negligible thickness along peripheral . direction. A 
single grain may have many tiny cutting points on the 
surface and during interaction with the workpiece may 
cause material removal. The transverse shape of the groove 
(Figure 1.10) can, however, be approximated by a circular 
arc jj+8j • Further, in most of the earlier studies 
{57 - 6oJ , it was assumed that there is no overlapping of 
cuts. For continuous feed, successive tips lying in 
different radial levels may produce cut at a particular 
point whenever there is grain-workpiece interaction. The 
generated profile height at that point will be the lowest 
height left after removal of material by cutting grains. 
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1 .9 Present Model 

In the present model , the grinding process is 
simulated to study the surface generated under different 
cutting conditions. The dynamic influence of table speed, 
and wheel depth of cut on grain distribution is analysed 
through the radial distribution parameter. The probabi- 
lity density functions of position of grain tips on the 
wheel surface are obtained on the basis of interference 
of grooves produced in transverse direction and the grain 
elastic deflection in radial direction. The regression 
equation of radial distribution parameter with depth of 
cut and wheel speed is obtained by comparing the power 
spectra of the experimental profile with the theoretical 
profile generated using Monte-Carlo simulation. The study 
is confined to single pass surface grinding under plunge- 
cut conditions. The transverse shape of the grooves 
produced by individual grains are assumed to be circular 
in section without side pile-up. Further, the grains are 
assumed to move in a circular arc which is approximately 
true since the wheel speed is much higher than the table 
speed. The wear of the grains are assumed to be negligible 
for a single pass. Only the grains distributed randomly 
along the radial direction within the effective profile 
depth, d , from the outermost wheel surface are considered 
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to be contributing to the surface generation and. not those 
lying beyond d . Considering successive cuts by grains 
along the cutting direction s longitudinal surface profile 
can be obtained from the grain- workpiece interference 
pattern. Transverse surface roughness and generated 
profile can be obtained theoretically by considering the 
interaction of grooves produced at different points where 
the groove mid-points are separated by a random distance 
called axial pitch. 

Once the radial distribution parameter under 
various cutting conditions is known , an approximate ana- 
lysis can afford a simpler way to obtain the undeformed 
chip thickness and the surface roughness. 
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CHAPTER II 


PROBABILISTIC MODEL AND SIMULATION 

2.1 Introduction 

The cutting surface of the grinding wheel can 
be thought of as a number of infinitesimally small 
cutting points distributed randomly over the surface. If 
the material is removed from a smooth surface, the result- 
ing surface will have a random profile characterising the 
grinding wheel surface. If the probability distribution 
function of the random cutting points on the wheel surface 
is mathematically known, the output surface profile can 
be generated knowing the input surface profile and input 
grinding conditions. Thus, the process can be modelled 
with grinding wheel surface characteristics as the black 
box having inputs such as grinding conditions and input 
surface profile of the workpiece and output as simulated 
profile of the surface. 

2 • 2 The; Axial Pitch 

The developed surface of the grinding wheel will 
appear as shown in Figure 2.1 a, where the various columns 
have been shown to be passing through the grain tips . The 
distance between the coloumns, defined as axial pitch, p, 
is a random variable. The tip position of ith grain in 
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jth column is shown in Figure 2.1 b, where y^ is the 

ti 

distance of ith grain tip measured along the periphery 

;t 

and Zj us due distance measured radially inward from the 
circle passing through the outermost grain tip. For any 
simulated width, , of the grinding wheel, the number of 
columns, H c , will be approximately given by 



where p is the average value of the axial pitch. Let L 
be the maximum width of cut produced by grain during 
interference with workpiece and fp (p.)) be the probability 
density function of the random variable p. The influen- 
tial width of groove is produced by any grain during 
interference and therefore, the pitch, p, will lie between 
0 and L. Hence, the average value of pitch, p, can be 
evaluated from 


P 



L 

j P f p (]?) ) dp-. 
»*** 


( 2 . 2 ) 


2.3 Probability Density .Functi o n 
2.3*1 Axial Pitch 


In Figure 2.2 a, (x, y, z) is a right handed 

0 

coordinate system fixed at the leading edge of the work- 
piece • At any section BB along the transverse direction, 
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grooves produced by any grain on 3 and 3 +lth columns 
(Figure 2.2b). As the grains pass independently of each 
other through the workpiece, x^ and x ^ +1 are independent 
random variables. The population under consideration for 
the independent section 0, L along axial direction has 
a distribution function F^ (x) and density function 
f-^ (x) (Figure 2.3 a) given by 


f X ^ 


I 

L 

0 


(0 4x < L) 

(x <C O 5 x />■ L) 


(x) =0 (x ( 0) 


x 

= ^ 

= 1 


(0 4 x 4 L) 
(x > L) 


(2.3) 


(2. If) 


Assuming the surface generated by successive interference 

of grooves produced, p, to be given by 


p = g j+1 - 

Using the transformation 


(2.5) 


w = x 3+1 " *3 

the probability density function of w is obtained 

C 62 l as 

f w (w) = / f x (»-*,> f X. ( - x 1> d x 1 

0 ■ - 


(2.7) 
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which is of t ri angular f onn (Figure 2.3 b ) . As the random 
variable Xj + .j °hd are uniformly distributed between 
0 and L, and f yj (w) is an even function 

fp (p) = 2 f lfI (w) U (p) (2.8) 

as shown in Figure 2.3 c, where U (p) is a random variable 

uniformly distributed between 0 and 1. Using equations 
(2.7) and (2.8), we get 

fp (p) = -§“ (L - p) (2.9) 

L 

Therefore, from equation (2.2), we get 

P = ^ L (2.10) 

Assuming the groove shape to be a circular arc of radius P 

(Figure 2.2), the maximum width of cut, L, will be given 

by 

L = 2y/(2f d 0 - d^) (2.11) 

where d is the effective profile depth. Neglecting 
the higher order terms in the above equation, 

L -r* 2 v/2""f"d e ' (2.12) 

For grain size greater than 36, it has been shown | 28^j 
that grain tip radius can be approximated by 

•p or g/2 (2.13) 

where g is the nominal grain diameter and is approxi- 
mated by 



34 


S 


0.6 

S 

g 


(inch) 


(2.14) 


where S represents the grain size. 

O 


Therefore, 
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3 Ve 4 e 


and 


3 
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\/ g d. 


(2.15) 


( 2 . 16 ) 


2.3.2 Grain Positions 

In any column passing through the tips of 
grains, the position of grain tip will he random because 
of the random occurremce of grain on the wheel surface. 

It has been shewn 54, 57 ! that the random variable Ay, 
representing the interval of successive appearance of 
grains for a fixed point on the workpiece, has an expo- 
nential probability density function, fy (Ay), given by 

1 - 

% ^ ~ ~J{ Q 


(2.17) 


wiie re, 


and 


AY 


A 


*5 - 


TTD 

M 


.( 2 . 18 ) 


(2.19) 


In equation (2.19), M is the total number of grains in 
any column. Therefore, M. is obtained by dividing the 
total number of grains existing on the wheel surface by 
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the number of columns, N . Hence, M can be approximated 




- T ' "D C 


( 2 . 20 ) 


■where C is the nominal grain density. From equations 
( 2 . 19 ) and (2.20), we get 


G/g d 


( 2 . 21 ) 


The radial positions of grain tips ( z ) are 
random. From the nature of the curve shown in Figures 
1 .4 and 1.6, the radial position of grain tips can be 
obtained f 29 , 33 j. Different curves are obtained for 
different wheels. These curves are generally extrapolated 
in the finer region where the values of t are very small. 
Various probability density function, f„( z ), of the 
random variable, 2 , have been used in the past jj>3, 54, 

58, 61 . From the nature of the curve shown in Figure 1.4, 

the general form of the probability density function appears 
to be of the form 


f z (z) = A z k (2.22) 

■where A and k are constants depending upon the cutting 
conditions and grinding wheel. Most of the forms of 
f z (z), used earlier jj?4, 58, 61 J can be obtained from 
the above equation assuming particular values of A and k 
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Under dynamic conditions, the outermost grain 
gets elastically deflected inwards f~61 | . The radial 
distribution of grain tip positions is thus modified. 

It is the conditional distribution of grain tip posi- 
tions In the radial direction (after the grains get elas- 
tically doilected) which defines the radial positions of 
grain tins on the wheel surface. The probability density 
function fl (z) and the corresponding distribution func- 
tion i'g (z) arc shown in Figure 2.4 a and 2.4 b, respec- 
tively. If d Q is the amount of deflection undergone by 
lowermost grain (Figure 2.4 c), then the grains will be 
lying between d Q and d Q + t along the radial direction 
of wheel surface. 

In Figure 2.5, ith grain lies at an average 
linear distance, ]\ » measured along the wheel periphery 
behind (i - 1 )th grain and the radial distance between 
the two concentric circles passing through the tips of 
ith and (i - 1 )th grain is Az. The.- maximum depth of 
cut by ith grain, is 

t ± = 2 A $ s/I - Ai C2.23) 

= t - Z\z (2.24) 

The condition for the ith grain to be an active cutting 
point (at least to contact the workpiece) is that the 
numerical value of t. should be positive j~63jj , i.e. ,Az 
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must lie between 0 and t. It has been found j 33., 61 |, 
that the value of d 0 is approximately equal to the 
maximum value of t in the fine grinding range. There- 
fore s for any grains lying between d and d + d , the 
probability density function of radial position of grain 
tips will be given by 


g z U) = f z a ! d 0 <i a D H- d e ) 

f 7 (I) 

= ' d +d ~~ 

,o . e 

f z (z) d z 
d 

o 

which on simplification becomes 


( 2 . 25 ) 


s z (2) 


(k + 1) (z) k 

(d + d ) k+1 - d k+1 
o e o 


( 2 . 26 ) 


In equation (2,22) and (2.25), the parameter k is a 
function of the cutting conditions . This has been gene- 
rally defined as the radial distribution parameter. The 
magnitude of elastic grain deflection, d Q , is generally 
of the same order as d g jj5^] . Hence, d c can be taken 
to be approximately equal to d g . Now using the trans- 
formation (Figure 2. 4- a) 

1 ii 

z • = Z’ + d 0 
8 


(2.27) 
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we get 


f , 7 (z ) 


(k + 1) (z + d fi ) k 

7k+i fjk+i _ ' ~ 
e ' * : 




( 2 . 28 ) 


Ttie probability density functions of random positions of 


grain tips on wheel surface in any column are expressed 


by equations (2.1 7) and (2.28). 

2 .4 Simul ation jof Surface Profile 


2.4,1 Transformation of Grain Tips to Workpiece 

The theoretical surface profile can now be 
generated using computer simulation techniques. The 
coordinates of the random positions of grain tips can he 
generated using the probability density functions in equa- 
tions (2.1?) and (2.28), 'The surface is produced due to 
individual interaction of grains and workpiece. Since 
the speed ratio is high, the path traced by an individual 
grain will be assumed to be circular as discussed in 
Chapter X. It will be further assumed that the material 
is removed either in the form of a chip or through plastic 
deformations without side pile-up, whenever grain-workpiece 
interference occurs. Depending upon the interference, the 
grain tip positions can be transformed to the workpiece 
co-ordinate. The locus of the transformed points on the 
workpiece during grinding will be. a circular arc and hence 
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the tip point of the grain after being transformed to 
Lli© workpiece will be the lowest noint of the circular 
arc « 1 : i e . j o kinematic translormations can be represented 

by 

y 1 = yVq - (2.29’.a) 

and 

z 1 = z" 1 (2.29 -b) 

The locus thus formed will be a circular arc of radius 

i 

R , where 

R 1 = R - z 1 (2.30) 

1 n 

having its centre at (y /q, R); R being the radius of 
the grinding wheel. Hie schematic arrangement is shown 
in Figure 2.6. From equations (2.17) and (2.29 a), the 
probability density function, fy ( Ay), of the interval, 
Ay, of the lowest positions of tips in the groove along 

the direction of cutting will be 

A i -At/A- 

fy (Ay) = -fq e (2.31) 

wiie re 

y\,j ss A-/q ■ (2.32) 

Using equations (2,28) and (2.29 b), the probability 
density function of the lowest position of ith grain 
tip in the groove along the vertical direction will 
be 

f 2 (z i ) SB fry (z ui ) _ (2.33) 
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Hence the lowest positions o.f tips in the groove with 
respect to the (x, y, z) coordinate system are 


i i+1 

X. = XT 

;i j 


M 




X^ 1 + L (1 -v/1 - R u ) 


i i-1 

yT = yt 


A i log (s^ ) 

1-1 


(2.34) 

(2.35) 


and. 


r~ 


d e | (2 k+1 - 1) F^_ + 1 


— i 1/(k+1 ) 


( 2 . 36 ) 


where. 


i = 1, 


M 


.1 = 1, . .., H 


and R , R and R denote the random numbers uniformly 
u i-1 U :L u 5 

distributed between 0 and 1. Hie transformations of 
probability density functions to random variates are given 

in Appendix A. 


2.4.2 Longitudinal Profile 


Once the lowest position of grain tips in the 
groove and its locus is obtained, the longitudinal sur- 
face profile generated by grains in any column can be 
obtained in the following manner; 

(a) The path traced by the preceeding grain may be 
overlapped by the succeeding grain as shown in 
Figure 2.7 a. In this case, the surface produced 



will correspond, to the path traced, by the succeeding 
grain and the active^ grain responsible for surface • 
generation is the succeeding grain. 

(b) The path traced by ith, i + 1 th and i + 2 th grains 
are shown in Figure 2.7 b. It may happen that the 
path traced by i + 1 th grain is overlapped by patii 
traced by i + 2 th grain. Though in this case, all 
those grains are removing material, actual grains 
responsible for surface generation will be ith grain 
and i + 2 th grain. 

(c) It may happen that the paths traced by ith, i + 1 th 
and i + 2 th grains arc such that none is overlapped 
completely (Figure 2.7 c) . In this case , all these 
grains will remove material and will be responsible 

for surface generation. 

This complex behaviour of the grooves produced 
in longitudinal direction by different grains is due to 
the random grain-workpiece interference . As a result, 
the number of active grains removing mate rial or producing 
chip and the actual number of grains responsible for sur- 
face generation under any cutting condition will be 
random. Combining the cases a, b and c, a general algo- 
rithm to generate the longitudinal surface profile can be 
developed in terms of the coordinates of the intersection 
points of the paths produced by the grains. 



k-2 

II I> (i, 3) represents the intersection of the 
paths braced by ith and jth grain in any column, then 
height of the intersection point z (i 5 3) will be 

given by 


( i 3 J ) 


1 

z + 


2 (R - z 1 ) 


( 2 . 37 ) 


where 


y 3) = y 1 + 0.7 


Xz lr.z 1 ) 2 B ( i _ I, 

t— (y - y ) _ J 

(2.38) 


Here y (i, j) and z (i, 3) are the coordinates of the 
intersection point B (i, 3) with respect to (x, y, z) 
coordinate system. The derivations of the above equations 


are given in Appendix B. Prom Figure 2 . 7 , it is observed 

that 


for case (a) : z (i, i + 1 ) ^> d , 

for case (b) : d z (i, i + 1) / z (i, i + 2), 

,y 

and 

for case (c) : d e / 2 (i, i + 2)^> z (i, i + 1). 

Extending the above algorithm for M number of 
grains in N number of columns, the longitudinal surface 
profile generated by the grinding wheel in different 
columns can be obtained. The computer flow chart for 
numerical computation of longitudinal surface profile 
height is given in Chapter IV. The number of active 



grains responsible for surface generation in each, column 
can 'also be obtained using the same algorithm. 

2J+.3 transverse Profile 

Once the longitudinal profile along different 
columns are generated, the vertical heights of the sur- 
face profile at different column positions along the 
transverse direction can now be obtained. These vertical 
heights will correspond to the mid points of the circular 
shaped grooves. Superimposing the groove shape by a 
circular arc of grain tip radius, /* , and considering the 
intersections of different grooves, the transverse pro- 
file can be generated. There are different possibilities 
of intersections of grooves. The groove with its mid- 
point at j + 1 th column may be overlapped completely by 
groove with its mid-point at jth column, (Figure 2.8 a), 
or is overlapped by grooves with its mid-points at j th 
and ,1 + P th column, (Figure 2.8 b), or none of the 
grooves are completely overlapped and have intersections 
with each other (Figure 2.8 c). The transverse profile, 
thus generated in each case , is shown by hatched lines _ 
in Figures 2.8 a, 2.8 b and 2.8 c. 

At any point p (x Q , y Q ) where x Q is the axial 
position along transverse direction and y Q is the longi- 
tudinal position along cutting direction with respect to 
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the fixed coordinate system (x, y, z), the resultant 
height of the surface profile, z (x Q ), can he obtained 

from the following algorithms 


z (xj = Kin j (y 0 > + <*o - i 5 


O 


k = 1 , •••? j 


(2.39) 


whore, 




Kin j d e , z£ (y Q )j i 


"1 , « » » j h , 


(2.40) 


and 


2 


, (y D - y k } 

_1 / ^ _ 2 - JL-~ + z 

Z ,, \Yr.^ ~ 1\ k 

k ° 2 (R - z k ) 


(2.41 ) 


In order to calculate the number of columns 
responsible for surface generation, a similar procedure 
as for bM,o mining the number of grains responsible for 
generating the longitudinal profile, is applied.. First, 
the intersection point of ith groove with rth groove 
is obtained whose coordinates are given by 


x 


1, r 


V 2 \ 


fOt_ (y Q > - z iflal 


1 ! '" f '"o' i ' J Qii g - + (x„ - x n . ) 


K| r __ "I ■■■ 1 1 ■ 1 

er r - v 


*i.r 


r - *4. 

z. (y D ) + 


r i 

(2.42) 

(2.43) 


g 


.here z ± (y Q ) denotes the longitudinal profile height 

in ith column at 


the transverse section at a distance y Q 
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from the reference coordinate system (x, y, z). Here. 
x i r-’ z i r* arG coordinates of the intersection points 
of ith and rth grooves with respect to fixed coordinate 
system (x, y, z). If z ± k is the minimum of all the 
Zi r , r varying from i + 1 to N c , then the transverse 
profiles './ill be produced by the intersection of ith and 
kth grooves. 

At any transverse section at a distance y Q 
from the reference coordinate axes, the increment Ai 
in the axial direction, is given by 

Ax = x m - x m-1 ; m = 1, N g (2.44) 

where H is the total number of sample points . The 

w 

transverse profile height at N g number of sample points, 
can be obtained using equations (2.39) 5 (2.40), (2.4l ), 
and (2.43). The transverse surface profile can thus be 
generated by joining these points. Computer flow chart 
for generating the transverse profile height at different 
sample points is given in Chapter IV. 
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CHAPTER III 


HEGKE-3SI0N AND SPECTRAL ANALYSIS 

3.1 Introduction 

In the previous chapter, it has "been shown that 
the statistics oX. the surface profile can "be theoretically 
determined by knowing the values of radial distribution 
parameter, k, and effective profile depth, d . Since these 
parameters are dependent on variables such as table speed, 
v and wheel depth of cut, d, the functional relationship 
can be obtained using the well established, methods of 
Design of Experiment and Least Square Estimates. 

Design of Experiment is a powerful tool for 

obtaining the functional estimate of the output parameters 

such ns k and d which can be represented as a function 

0 

of the input parameters such as v and d. The least 
square estimates of the output parameter can be evaluated 
at the minimum of the sum of squares of the difference of 
trie power spectral density of the simulated and experimen- 
tal profile. The workpiece profile is the sample realisa- 
tion of the stochastic grinding process and its statistic 
is characterised by its spectral density or autocorrelation 
function. The effect of v and d on radial distribution 
parameter, k, and d can be estimated using analysis of 
variance. Ap flying regression analysis, they can be eva- 
luated as a function of v and d. 



3.2 Design of Experiment 
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3,2.1 Fac to r i al Exp e rim e nt 

Design ol Experiment is defined as the procedure 
for selecting the number of trials and relevant conditions 
essential and sufficient for solving the problem that has 
been set at the required precision. Hie output par am eter, 
k, can be mathematically represented as 

k = 0 (v, d) (3.1) 

■where 0 is the response function and v and d are input 
parameters known as factors. Each factor can take several 
values , called levels, in any trial which can be obtained 
by treatment combination of factors. The selection of the 
experimental region of factor space is connected with a 
priori information. The variations of v and d in fine 
grinding are usually in the ranges of 5 to 15 m/min, and 
4 to 12 yU (micron), respectively. The response surface 
of k or d is shown in Figure 3.1 • Ihe region of deter- 
mination of the factors v and d can be estimated by 
taking a section of the response surface parallel to the 
piano (d o v)and projecting lines, thus obtained on to 
that plane (Figure 3.2). An experiment in which all the 
possible combinations of the factor levels are realised is 
called a facte; riel experiment. As the factors, in our 
study are two, wo have a 2 x 2 factorial experiment. 
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3.2.2 Repetition of Trials 


The running of repetitions of trials never gives 
coincident results since some experimental error (error of 
reproducibility) is always present. This error is esti- 
mated by repeating the trials. The variance of k in i*th 

p / \ 

trial , Si, (It) consisting of repeated observations is 


calculated using the formula 

1 


3?, (k) 


n 


1 


_ 2 

U=1~" 


(3.2) 


for i* = 1, ..., N-j . Hence i* denotes the trial having 

treatment combination of ith level of one factor and 3 th 

level of other factor, N 1 is the total number of such 

trials, and is the ne an value of k obtained in 

repeated trails given by 

n. 


*1 


k 




X_ lr 

t: k iju 

U 1 U=1 J 


( 3 . 3 ) 


■where k 1 is the value observed at uth observation for 
the i*th trial. The variance of reproducibility , s rep 
of the parameter, k, can now be obtained from 



rep 


(k) 


~ lij (n 1 -17 



i’=1 


n 1 

.c 

u=1 




(3-^) 


for equality of the number of repetition at all the trials, 

otherwise, from 
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N, 


°rop Ck) 


N. 




f i> s 4 


( 3 . 5 ) 


i=1 


where f ±1 is the number of degrees of freedom in i'th 

t/ X'X Q. JL j 3- # 0 * j 


f i' = n i - 1 


( 3 . 6 ) 


3.2.3 Homogeneity of Variance 


One of the requirements of regression analysis 
is the h or, one unity of variance which means that there are 
no variances among the ones being summated that consider- 
ably exceeds the remaining one. As (k) and S . , (k) 

T9p 1 

are of the same order, homogeneity of variance can be 

assumed >65^. When it is thought that the population 

+•*! 

variances are not equal but more of less randomly scattered, 
the test most often used is the Bartlett's test j 65]. The 

test statistic 


N, 


jL 

c. 


? log ^rep 


(k) - f i' log S i» 

i =1 


(k) 

0 . 7 ) 


when 


c 


2 


0 . 434-31 1 





( 3 . 8 ) 



f 


and. 



i'=1 
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c 2 approximately obeys a chi-square distribution with 

(N,| - 1 ) decrees of freedom. For ensuring homogeneity 

of variances, the value ,X r f (N. - 1) obtained from sta- 

v '"c 1 

tistical table with (Ih - l) degrees of freedom for cL % 

* o 

significance level, must be less than B. 

3.2.X Analysis of Variance 

Through analysis of variance, relating to the 
sum of squares identity, several sources of variations 
can be isolated, estimated and tested. The least square 
estimates of k, for the randomised 2x2 factorial ex- 
periments have been given in Table 6.3. The values of d Q 
obtained from the experimental profile (Figure 3.3 a to d) 
arc listed in Table 6.1 . The analysis of variance table 
for k or can be prepared with the following identi- 
ties and their corresponding degrees of freedom |_66j . 

The identities and degrees of freedom are 

SSd » Sum of squares for factor, d 

i . y 1_ " T 2 _£ 

- ^ a, hi 1 " °1 r 1 n l 

SSI = Sum of squares for interaction 
S SS Tr •• - SSd - SSv 


(3.9) 

( 3 . 10 ) 
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SSv = Sum of squares for factor v 

r. 


1 


c 1 


1 


5=1 


2 

I , 

• .1 * 


T £ 


C 1 r 1 n 1 


( 3 . 11 ) 


sss 


Within sum of squares 


SSI - SS Tr 


( 3 . 12 ) 


S8T 


Total sum of squares 


* 

i=1 


\ 


5=1 


n l 



u=1 


i 2 

K iju 


C 1 r 1 n 1 


( 3 , 13 ) 


SS Tr = Treatments sum of squares 

- C 1 , r l 2 ’ 

1 5 > T. . - * 

i=1 3=1 


fr 


n 


1 


C*j P-j 


= Duqreo:: of freedom for d 

a c 1 - 1 

= Decrees of freedom for interaction 
= (c^ - 1 ) (r,j - 1) 

s Dt -:,u Hi a of freedom for within 


c 1 r 1 (^ 


1) 


Do i -trees of freedom for total 


c 1 r n n 1 


1 


( 3 , 1 ^) 


( 3 . 15 ) 


( 3 . 16 ) 


( 3 . 17 ) 


( 3 . 18 ) 


: Tr 


Decrees of freedom for treatment combination 


f d + f v + % 

Decrees of freedom for v 
r- - 1 


( 3 . 19 ) 


(3 


1 



58 


In the above 

enuation 

S 




n-j 




i j • 

"u~1 

k i t 1u 


( 3 . 21 ) 


. r 1 

n 1 



T. = 

I# * 

_ 

,1=1 

• 

... 

U=$ 

k. . 

13U 

( 3 . 22 ) 


0 *j 

Sl 



T = 

• j » 

i=i 

"\ 

Z_. 

u=1 

k. . 

( 3 . 23 ) 

and 

c i 

— 

T=f 

r 

4 1 - 

s' k_. . 

iju 


T 

y 

T=T 

( 3 . 24 ) 


Hie values of the above identities for d and k are 

G 

given in Tables 6. 2 and 6.4-, respectively. Prom Table 6.4-, 

it can be seen that the computed F - ratio for d to 
’•within*' is larger than the value obtained from statisti- 
cal tabic, for the upper *>% level of the value of F 
distribution with the f^ and f g . The statistic F - ratio 

for d is given by 

F - ratio = IP / ( 3 . 25 ) 

1 d Z Q 

Thus, it uan bo concluded. that the effect of depth of cut 
on the j-n.ci.Ln3. distribution parameter, k, is non-zero. 
Similarly, the effect of table speed on k can be obtained 
and which is also found to be non-zero, but the effect of 
interaction of d and v on k is negligibly small §s- ' 
the- Co* spa tod F ratio is very much less than the value 
obtained from statistical table. Proceeding in the same 
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manner, the sum of square identities can be obtained and 

hence the analysis of variance table for d (Table 6 . 2 ) 

6 

can be prepared. It can be seen that all the computed 
values of i 1 '- ratio are much smaller than the theoretical 
values obtained from the statistical table. The effect 
of d, v, and their interaction on the effective profile 
depth, d , can therefore be assumed to be negligible and 
hence d can be taken as the mean of all the values 
obtained expo rimen tally . 


3 • 2 . 5 Least Square Estimation of k. 


The least square estimates of k can be obtained 
from the factorial o'-yortnont at the minimum sum of squares, 

r~ ~~j 

SSk, fay Riving iterative values of k } 62 j from the follow- 
ing equations 


m ^-1 


SSk 


Min 



(S' (fq,T) -S (fq,S- , k» (3.26) 


h=° 


^1 fC 

at different frequency, fq, equal to — >Aiere it is 


“l 


assumed that 


c" (fq, "if ) 


'S (fq, % , + 


f(fq) (3-27) 

In the above equations, €(fq) is uncorrelate* erro r 
at frequency, fq, 1 represents kinematic and operating 
parameters (v, d and Sg), G (fq, 5 ) is tlle smoothed 
spectral density of the actual profile at the frequency, 

fq.for Mrindinn condition, 1= , ^ (fq, ^ > k) is 
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smoothed spectral estimates of the theoretical profile 
simulated at; the grinding condition, "> , assuming some 

j 

known values ol k and the number of spectral estimates 
over a frof money range of 0 to cut-off frequency, fc, 
is taken urmnl to the maximum lag number. 

3*3 Syoc t .ml Arialv s o s 
3.3,1 Power Spectral Density 

A single time or space history representing a 
random phenomenon is called a sample function (or a sample 
record when observed over a finite time or space interval). 
The collection of all possible sample functions which the 
random phmmm'.uioi ; might have produced is called a random 
process or stochastic process. Hence a sample record of 
data may be thought of as one physical realisation of a 
random process. The workpiece profile produced by the 
grinding wheel is a sample realisation of the random grind- 
ing process. Power spectral density function measurements 
of physical data establish the frequency composition of 
the data which, in turn, bears important relationship to 
the basic characteristics of the physical system involved. 
The mean square value of a sample time (or space ) history 
record in a fro (money range between fq and fq + Afq may 
be obtained by filtering the sample record, z (1), with 
a band pass filter having sharp cut-off characteristics, 
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and computing the average of the squared output from the 
filter. Thin average squared value -will ...approach an 
exact mo an square value as the observation time, T‘, or 
length, L, approaches infinity j 67J - Mathematically, 

this can be expressed as 

2 ' / L 

Y [ illj i’q + A fq j = Lim J / z 2 (1, fq, /\fq) dl 
f z - L L A ’ 


( 3 . 28 ) 


where z (1, fq, fq) is the portion of z (l) in the 

2 f— — 

frequency range fq and fq. +/\fq and A [ fq, fq + Afq 

' z i_ „ 

is the mean sonar c value of the power spectral density 
from fq to fq + A fq. 


For small A^q, a power spectral density, 

G„ (fq), can be defined such that 
z 


r\ [fq. fq + Afq] ~ o 2 (fq) Afq 


( 3 . 29 ) 


Therefore, 


G z (fq) 


Mm ^ z fq + A f H 

£fq-+0 A f< l 


Lim Mm 1 — ■ / z^ (1* fq,Afq)dl 

Afq-^O L-*c* WL y Q * 

( 3 . 30 ) 


For stationary process, the power spectral den- 
sity is related to the autocorrelation junction, R z ( <T)j 
of the process by a. Fourier transformation as follows 
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ir-0 


, r .. -.1 2irfqr 

' : \ C - ) G dr 


(. (i'M) 


R ( • 


\ . 
*}' 


/ il.. Cl) cos 2Tf fq T dT 


(3.31) 


• mal valued, non-negative function as 
n.v! even function; 


3.3.2 Aulooorro l -lion p'u notion 

'•‘I*'-- e u toco rroln Hon function for random data 
cl o sc ri i» ts! e* g moral dependence of the values of the 
(into ;*!, no i a I, on the values at other point. An esti- 

mate for the ••n. (^correlation of a (1) at any length 1 
and 1 + 7 : ny 1 ift obtained hy taking the product of the two 
vnlui.-:: nee .••vu-r.nf, i a/- over the observation length L. The 
r* us ilk i eg c,v<-rnp.- nroduet will approach an exact auto- 
corre in Lion pi-oiinrt ns L approaches infinity. The auto- 
curr* i n !..i ot: function, defined as R ( L ), is given by 


K w ( 7 5 l ..v, 


J- f z (1) z (1 +T) dl (3.32) 
t) 


(7’) is always a real valued even function with maximum 
iv bo cither negative or positive, i.e., 


at [ ~ 0 and m. 


1 (70 = 

and 

R ( * ; ) > 
z yy 

for ;*i l C , 


R„ (-71 


ii, 


(3.33) 


(3.34-) 
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The mean value of the process, can 6e expressed in 

terms of autocorrelation functions as 

= VR Z (jg)' (3.35) 

.2 

Similarly, the mean square value, r ,can he given by 

\p 2 = R z (0) (3.36) 

i.e., mean square value is equal to autocorrelation at 
zero displacement lag. 

3.4 Multinle Linear Regression Analysi s 

If the criteria for homogeneity of variances at 
different observation in any trial is satisfied and the 
levels of the factors are non-random and uncorrelated, the 
multiple linear regression analysis can be applied to esti- 
mate the regression coefficient. In this analysis, the 
random response of k is assumed to be linearly dependent 
on d, and v because the effect of the factor d and v are 
different from zero as discussed under analysis of variances 
(T a ble 6.4). It is further assumed that each observation, 
k^, in any trial, is randomly drawn from a normal population 
with mean, • , where denotes the true regression of 

on d and v and given by 

n ± = 1 0.37) 

for i = 1, ..., N 2 and variance ^- 2 . A> A 


are 
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variable parameters, and v^ are wheel depth of cut and 
table speed, respectively for ith observation and is 
the total number of observations. Estimators a^, b^ and 
of , and P ^ are found so as to minimise the 

square identity, Q, given by 


Q 



Z 

i=1 



b i a ± 




( 3 . 38 ) 


Setting the partial derivatives of Q with respect to a 1? 


and b^ 

equal to zero, 

we get the normal 

equations 



3. 

V— 2 


H 2 



1^2 

+ ( d, ) b 1 

i=1 

+ ( Z T ± ) 

l=l 

b 2 = 

It 

1=1 

( 3 . 39 ) 


N 2 

IZ 

Z 

) 

h 


M 2 


( Y 

1=1 

d. ) a. + ( 5 

1 1 1=1 

d?) b 1 + ( 

1 d. 
i=1 

L V i^ 1 °2 

= > a 
1=1 1 







( 3 . 4 - 0 ) 


Jz 

n 2 


N 2 


n 2 

sr** 


(2_ 

i=1 

v.)a 1 + () > 

1 1 1=1 

a i T l> b 1 + 

( > 
i=1 

v i > b 2 

= Z\ v i 

i=1 

^i, 






( 3 . 4 - 1 ) 



The solution of the above equations gives the unbiased 
estimators of the parameter ft ^ and ® le estimated 




regression, of k^ on d- and v. will be 


k ± 


a i + b i a ± + b 2 v i 


( 3 . 4 - 2 ) 


2 ' 

The unbiased estimators of (f can be shown to be 


N, 


S; 


k/d,v 


N, 


2 _ ^ 

(k, - Y) 


( 3 . 4 - 3 ) 


1-1 



whe re , 

<<> 

k 


JL 

N, 


* 

i=1 


k i 


Using 


d = 


v 


and 


k 


JL 

lh 


JL 

N, 


JL 

N, 


d i> 

1=1 




V i’ 


*2 


i4=1 


i» 
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(3.44) 

(3.45 a) 

(3.45 4) 

(3.45 c) 


the normal equations C 3 • 39 ) » (3.40) and (3*41) can be 
further simplified to the form 


a 11 

b 

+ a 12 

b 2 


a 21 

b 

+ a 22 

b 2 

= s 2 » 

and 





cl* J 

SS 

5 - d 

*1 - 

’ v b 2 

idle re 





a 11 

=S 

fa 

> 

i=1 

(d ± 

C\] ^ 

l 1 ^ 

1 



N 2 


_ 2 

a 22 


H 

(v i 

- v) , 



i=1 


No 

a 12 

= 

a 21 

- 

Z (a i 





i=1 



N 2 



s 1 


£ 

i=1 

^ d i 

- d) (k. 


(3.46) 

(3.47) 

(3.48) 

( 3 . 49 ) 

( 3 . 50 ) 

- d) (v ± - v), 0-51 ) 

- k) C3.52) 
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and 


N, 


i=1 


(v± - v) (k ± - k) 


(3.53) 


Solving the equations ( 3.4-6 ) 9 (3. 1 +7) and (3 -4-8 ) , the 
values of a-j , b^ and t> 2 can be obtained and the regression 
of k^ and and v^ can be obtained from equation (3.*+2) 

after substituting the values of a 1 , b^ and bg. 

o 

The unbiased estimator of $ in equation (3 .^3) 

can be simplified as 2 



i=1 


SSk 


°k/d,v 

tiihere, 

SSk = 


t . h • 
3 i l 


N, 


N, 


~J=i 


3 


( ) 

i, ■ , — 1 1 . ... 

KL 


(3.5^) 


(3.55) 


Since k-'s are independent, the variance of b is 
i 

can be expressed as 


J b. 5 

i 


% 


2 


J n cr 


(3.56) 


where c_., can be obtained from the matrix !c_. 


'ii 


L ID 


3 


which 


is the inverse of the coefficient matrix a^j 5 obtained 

from equations (3.^+6) and (3*^+7). Using the unbiased 

estimator of 2 given by equations^ (3. 55 ) and (3.56), it 

can be shown that the statistic ( “§ “-Z ) has student-t 

Vdjv^ii 

distribution with (Ng - 3) degrees of freedom. The sta- 
tistic can be used to test any hypothesis of the form 



67 


= constant. The (100 - a - ) percent symmetric 

j -I O 

p 

confidence interval of is given by 



V 


/ 2 


O'L - 3) 


S k/d,v v / 




/ # 
<r± 



t ...r / 0 
- A c /Z 



3) 


S k/ d, v 


\/ 



(3.57) 


for i,= 1, 2, where can be obtained from the confidence 

interval. From equation (3.57), 95% confidence interval 
for fe-\ 3 -r.d can be obtained. 
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;• . ar !v 


IV 


• •• • i CAL b( if !TIOM TECHNIQUE 
4- . 1 I ntrofh u ~ 4:1 > >?: 

iv, oi' random data such as the surface 
height of vo v-:; «<>•;.. U e requires digital technique. 

Discrete rcr-rv::i ;ut-, i, i no oi’ continuous data involve problems 
of dieiti/C't i i‘0‘ l-iconfvtl; of integrals with sum and rela- 
ting record i on; .*> frequency band widths to correspond- 
ing (.!i: if:'. I :r r,- r. Ail these techniques must ensure 

the ."least 'ions <..f inform-sloti. The experimental surface 
profile is •; aiivus revs’ll uhj.cli needs sampling and 

digital ■ A:: . 

<:■ « •>•?!».• j* A coi-.-ru-o tfic experimental results with 
the them re ;.ir-.i •.*;<■- !. , simulated profile must be generated 
from U;v -• ;•< i-iiltv densities and statistical para- 
mo to.i*! iM! .f * * i * :* I tcciiniquos. 


4- . 2 111 1 i . t i ; 

continuous f ■ 
main parts, c 
The snj;ipli.;i, 
the are 

too far r* 1 1; ,..j 'i 


■Lh.is o i Continuous Data 

7 i , roc a. : ■ s of digitizing consists of converting 
i. discrete numbers. These involve two 
.; ■(> v- *•;;>, : r tL.ing and the other is quantization. 

•; -• n of i. ned as selection of the points at which 
t ; ; rv < (i . 1 f th (•; sampled value s are sepe rated 
1 >]. c.'d values would represent either too 
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low or too high frequencies in the original data. This 
property is called aliasing j6 7,] . If the space interval 
is (\ 1 between samples, then the sampling rate is 
samples per interval. The useful data will be from 0 to 

A 

2 7 ^ 2 cycles/interval (c/interval) , since the frequency 

-j 

in the data v/hich are higher than c/interval will be 

folded into the lower frequency range from 0 to and 

. confused with data in lower frequency range. The cut-off 


frequency 

fC = (4- .1 ) 

is known as Myquist frequency | 67 ]. One of the practical 
methods is to select fc to be few times greater than the 
maximum frequency of interest .which':, can be determined from 
the continuous record of the data. This. process will lead 
to choose Al small. If - A 1 is too small, then samp- 
ling will yield correlated and highly redundant data. For 
proper Quantization, it has been shown [67_j that the sample 
points should be greater than 256. Our values of sample 
points satisfy this criterion. 


4- . 3 Autocorr ela tion Functions and.. P pjJjLr.^^ 

4-. 3.1 Autocorrelation Functions 

The autocorrelation function can be obtained 
through the numerical calculation in Steps 1 to 6. 
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Sfco'P _1 " ^ '■ * ! • • '“to. * * 1 * selection of sample size. 

i;,; ' A t.! ( . , ) ' 4u; values of any sample record 


(1) 


, I 


>' 31 


them z O | . i i 

(-• ! J 

for n’ * 1 , ... 


! 


o 

id : s ; ro 


n ! A 1) 


(4.2) 


(4.3) 


h r g is the sample size such 


A 1 


and 1 in r ■ i aimer: 'oooro length. 


(4.4) 


■6ca n : * ? 

>■ M ’ A' J ' 

r; f 

i*k A 

, Sc.. i . 

r i i o 3 

la;-;, n ;j , in. 

'.-■I . •; < l.i 

! H : J.O 

a; hu*, Tii. 

i'uJV 1 

i.ho i.i.'J, 


the maximum number of correlation 


1 


A i 


(4.5) 


Ion "■> Cal mil oT the mean value. 


The r, amnio mean, z, is estimated from 


; i 

O 


1. z (1 n f) 
n’=1 


(4.6) 
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Sten 4 • r • r-i 


w zoxo mean value. 


u 




! i V. i / 


(1) -w.il.; 


v i . f / 


'■ ' X! iii story record u (1) where 

, m - 


>1 (1 ,) Given by 


•J. (:\ rt + r.' Al) 
( 1 ) ~ « 


(4.7) 


(4.8) 


Stop 


r " «. 
) *> 




if the U‘ -an square value, 


•r: square is Given by 




■> <V 


ft. 9) 


The; n s.k A ! I * i t;.' 

/yl V aa ; i ;.y 


A.rKJud r'jihima-hor of the true mean 


sample* standard deviation <T , is 


1/2 


I it, - 1 


X 

\ 

/ 


n=1 


u (1„) 
n | 


(4.10) 


Oton Q, OrU eu! - J.u.n of ; * u tocor.r elation function. 


if, data values u (1 .), n* — 1 3 from 

*1 n ^ 

the tranr: to used record, u (l), which is stationary with 
mean, ti =-. o {.he* autocorrelation function is estimated by 


the formula 


U 


it., (r ,Al) 
n 


11., - r 

i3 


Vf 

A u ^n’^ u ^n'+r^ 


n*=1 


(4.11) 


lor v 


ij 5 *3 * * * ? ll 1 
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If. 3. 2 Pover *>i tr-fl. Density 

:i • :l, ‘ tlo* i'Uy random process, the power spectral 
density ho. inf. the Fourier transform of the autocorrelation 
function oi' Fu: sample record, the raw estimate, ~(T (fq), 
of a true spectral density, G (fq), can he obtained from 
the dig! tal V’«'1 m-r; oi' autocorrelation function for an arbi- 
trary frequency fq in the range 0 g^fq fc j 67 ! from 

rn^-1 

/\ -t-T* — 

) : 


T (fq) = . 2Al \\ + 2 / R r cos ( II f~ :q 


+ H cos ( 
m l 


r=1 

tt m 2. 
fcT~ 


(4.12) 


The values of the .function G (fq) need to be calculated 

z 


only at the in.. + 1 spatial discrete frequencies where 

fc 

fq 


U 1 


m. 


(4.13) 


1 


for k 1 ~ 0, 1, 2, . m 1 . Ihe spectral estimates at., 
points loss than apart will be correlated. 


Step 7 : Haw estimate of power spectral density at 
harmonic k^. 

The raw estimate of the power spectral density 
function, G. , at harmonic, k-, , corresponding to frequency 

JVq 

k fc ^ i— — i 

fq = -1—-— can be obtained from \J>7 \ 
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i. c 


A 


ra-, -1 

J . 


li + p 
r, 


£ 

r=1 


Vv -IT r k-, 

H cos ( ) 

i m i 


+ (-1 


k-, /' 

n _L «n 


m 


i i 


(4-. Ilf) 


o!.' ! ‘ . ' 


ill,' o.f 


o.'ii.lni-i to of power spectral density. 


tiono 


rt-ci idod to chock the previous calcula- 
■ . ' • 1 i t ;i. o , 1 vo . i by j 67 | 


V 


i 


1 1 /rr" X > — i 

vr T i 5 L o + ^ 5 , \ + 2 \ j 


(4.15) 


.i f.’ 1 * i r*‘ v - J.. j 1 i ■ v oi' i nput data . 


Vi. : i ;-o*m of pro whitening, sometimes recommended 
i a s.'i . - •; f y ; ; i , amounts to applying a special filter 

(biyw.-j • i-,j if. ) ’,<■ Mio data which will result in the 

»VI ! r ' ■; t,.-.. Dab spectrum. This process can be 

aohiovfci ’.V If. -a? following moving linear combinations 

u (rrO = (z (l n , +1 ) - a) - c 4 (z (l n ,) - z) 

(if. 16) 

for n ,! j= 1, . (n 1 - 1), where c lf _ is a constant to be 

supplied and its value is 


* 

i 


i 




S i.o. 



(4.17) 
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While applying the prowhitening process (4-. 17 ) 3 Step 9 
is to 'b «4 onlculatad along with Step 3 and all the numeri- 
cnl SU.-nr; from l j to 8 with u (n’O in place of u Cl n ,) 

Step 1U ; Smoothing of the raw estimate of power spectral 

density. 


density 

spec tr* • 


Smoothing of the raw estimate of power spectral 
can be clone by Hamming [_68j. The smoothed power 
f i t >i i / f; N will be given by the following 

iOi 


enua 


i \ 

' ‘ O 





0 . 9 H G + 0.46 G ,j 


(4.18) 


0 . 23 ‘Cf _ n + 0-54 G k + 0.P.3 G k 


1 " 1 (4.19) 


for O K j y ia-j , • tiul 




o.H4 (< m _ + Gji-6 G m 
J. 


(4.20) 


•llm ri; i I il.y of t.lio estimates of the true spectral 
(k! .,. ;U ,y < ;.c, t ■;> .to!, t!ec reuse with increased record length. 
IVia I'i Stic can Ue achieved by smoothing. 

Stem 11 : Compensating or recolouring of the smoothed 
spectrum (if prewhitening is applied). 

Tile smoothed power spectrum, if prewhitening is 

applied, in compensated as 


G C iv^ ) 


- / (i" / (1 + c 2 - 2 ci*. cos (~\r y ) (lf,2l) 

k l 


for k-, = 0, 1? •••> rn i* 
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4.4 Long.H ' j. ■ : 1 nonsvMrso Sur f ace Profile 

bra 'Mur::; of h are estimated at different 
grin < ‘ i t i on-:, lie* longitudinal and transverse sur- 
face: • *r* *f.: U- in , ,m nor?,, tod theoretically through the 

numerical ov-. i ..-.f : n;i cm’ two following steps. 


Ston 1 » m* • • t; on c, f Li iron sets of random numbers uni- 

fan. .< ■/ distributed between 0 and 1. 

; 1 m j iv<r».:iy h i strlbutod random numbers lying between 
0 arici 1 t.,w > '< ; no v- 1 Loci by nneudo- random number generation 

method or v/i th the help of built-in- subroutine programmes 
HA rib v 1 (a), Iwl -Vrp(X), IiAMDI3(X), available in IBM 7044 
Computer. 

stop 2 ■; Of. i c,:ul'it;.ort of L, U Q , M, X, /\.j, and k 

• ■or ■*.!:, y r.pood ratio, q, grinding condition, , 
of foct.i. vc: n rofi In depth , d^, , obtained i rom experimental 
[JiuCij wrKi simulated grinding width, .>*(,? of grinding 

wheel of o i D, ! awing grain size, Sg, nominal grain 

density, C, the value of L, tJ,, A » A«» anci k are cal “ 
eulated fro- the.? equations (2*12), (P.20), (2.21;, (2.32) 

a nd (3.42) re s p e e tive 1 y . 

Stan 3 s Generation of random positions of lowest grain 
tips in any groove along cutting direction. 

Using the equations (2.34), (2.37), (2.36) and 
the random numbers generated in Step 1, the random posi- 
tions of lowest grain tips in any groove along cutting 

direction are generated . 
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Stop 


1 1 \j< h : l‘" 


i<y t \y. !. i *1 If Ill'll 


- , l 


m3, surface profile 

\'\'i >v. ) I. II iiii IS • 

. . ... rt .... ... no:;i tions of grain tip 3 at any 

. . , !, . ,iv itii column is calculated from 

soiU'Oo : , I ' 

the on ■•.••tic..-. 

i £ ~ 

( 4 . 22 ) 


•j. 

li 

i . VU&- 

= *]*;, (B-zb 

J 

where 3 

i "| . . .. I 1 1'/ I 1 ! i.'l.t C ’ 

‘a’ 


rioun sr-n/ersen in one. pass 

to vorh 

i f n 

.. profile height at 

f’i f » • ’ ' 1 



* ; () f d and all the values of 

Z, (l js ,) VJ.iU t:c tl-U ° 

' J 1 |v|. This is continued 

,/i ) nbi/ i .luoa ior i. - h ***’ 

'"j "’ n * . p and for all the 

pc . i 'i ror n* =1, — > 8 

Therefore the longitudinal pro- 

euli nuns, .1 - 1, •••» ‘o* 

, •, r . nP he obtained, 

fil.. height :.t din-.u-oni. eoiiii-'ii^ tin 

. . .. i. r-m-vorso profile at any section BB 

« lt . n K, Urm mtxori oj , 

r'-i Tection at a distance y 0 » 

, a .ong transverse collection 

• i 4-r* of all the columns 
Longitudinal profile height. 

. . •, from Step 4 substituting 

at a. uis sauce y Q groove with its 

, , \t The z coordinates ox any gr 

1 oo.ual to y c * * ' is 0 ht a ined 

. t - ( y ) n t the sample point n « 
midpoint at 21 j ^ o J 

from the equation 
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^ a , > (l n" ~ x i ^ 

z j (1 n” } = z j (y o } L (lf - 2 3) 

where , 

l n „ = 1 o" + n " A 1 " (4.24) 

and 1" is the increment equal to the simulated width 
dC j of the grinding wheel divided by sample size, N g , 1 „ 
is an arbitrary initial value. 'The transverse profile 
height at the sample point, z (1 will be the minimum 
of d and all the values of z. (1 ,,) obtained for j = 1, 

G o i-*- 

N * This iscontinued for all the sample points 
n n = 1 , . .., Ng and thus the transverse profile height can 
be obtained. 


Step 6 : Calculation of average number of active grains 

responsible for surface generation. 

After the Step 3, the coordinates of B (i, o ) in 
any i’th column is obtained from equations (2.37) and 
(2.38) for i = 1, and j = i + 1, to M. ' If the minimum of 
all the values z (i, j) is z (1, k), then the count of the 
number of grains responsible for surface generation upto 
y (1, k) is 2. Changing i = k, and j = i + 1 , 
the above procedure is continued and total count of the 
number of grains responsible for longitudinal surface gene- 
ration, (i’) } in i’th column can be obtained, following 
the same procedure for i* = 1 } y c > the average value 
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of the number of active grains responsible for surface 
generation in any column can bo obtained from the equation 

:; c 

— } h (4.25) 

! ^ ".J 


1 




c 1 

The longitudinal profile heights along different 
groove mid -positions at a distance y , are obtained after 
executing the Stop 4 replacing l n , by y For this trans- 
verse section lying at a distance y n , s the intersection of 
ith groove with jth groove is obtained using equations 
(2.42) and (2.43). This process is repeated for i = 1 and 
j = i + 1 to N . If the miniiiium of these values of . is 
m , then the count of the number of grooves responsible 
for transvors surface generation is 2. Changing i = m 
and j =111 + 1 to M and repeating the above process, the 
number of columns responsible for transverse surface gene- 
ration, N, (y , ) will be the cumulative sum of all the 
previous counts. Taking more.. number of sections, the 
average value of the active number of columns will be 


Wr 


N, 


4 


1 

N r 

J 


51 \ h n , ) 

n'=1 


(4.26) 


where. IV is the number of transverse sections. 

b 

Therefore the average value of the number of 
active grain density responsible for surface generation 


will bo 



8 ?. 


C. 


IT . \ 

rc Dc* 


(4.27) 


F nw Chart 

The computer flow charts for the above steps are 
shown in Figures 4.1 and 4.2. 
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5.1 Introduc tion 


CHAPTER ,V 


APPEOXE UIS AIIALYSIS 


It has been shown earlier that the radial distri- 
bution parameter 3 k 3 and effective profile depth, d , 
describes the distribution of grain on the wheel surface. 
In the simulated study of the grinding process d g was 
evaluated from the experimental profiles where it was 
assumed that the effective profile depth, d , is equal 
to the maximum peak to valley height of the work surface. 

When Interest lies in determining the values of surface 

/" 

roughness produced in fine grinding, a simplified ipddel 
can be proposed. 

— -Determination of Peal; t o Valiev He igh t- 

The surface produced on a stable grinding machine 
consists of scallops left on the workpiece due to finite 
spacing of successive grains. If all grains are assumed 
to be at the same radial distance and equally spaced, the 
plan view of the finished surface will consist of array 
of slender shaped scratches as shown in Figure 
Assuming the mid-point of the scratch to be the deepest 
point located at a distance h, defined as idealised mean 
peak to valley height , Jeri can be evaluated from 



h 


v 


(5»r) 



H 


',/ 2 /' D 


The value of C in the above eauation corresponds -to 

2. t 

z' = h (Figure ..1- .4- I and equation CJ.T) can be solved for 
h by iteration. As the grains are randomly distributed,, 
the assumption that the active grains extend the same 

distance from the wheel surface is a crude approximation. 

3 

Figure shows that the peal: to valley" height depends on 
the coarsest and finest spacings h^ and h respectively 
and upon the scallop height h g at z = h u L33^J* From 
this figure, it is evident that the grains producing 
grooves lie,:- within h u and generate the transverse pro- 
file. For any grain to be active, the grain depth of cut 

should be positive i.e., the grain must contact the vork- 
J3> 

piece jj>3"j. Under thiso- conditions, the limiting value of 
the radial zone within which grains contribute cutting 
is approximately equal to 


2 . 

h u = 2 7 VS 

where is the average linear grain spacing. Therefore^, 

the maximum pe ak to valley height of the surface can be 

taken to be equal to the maximum value of the undeformed 

/ 

chip thickness. If C-y is the active grain density then 
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/ 



1 



<5-3)’ 


where b is the average width of the groove produced by 

_J2_ 

a grain and can be obtained from \ 2 &\. 


I* 


b £ 


8 d 



( 5.^0 


% 

*/ 


Therefore, from equations (5.2) and (5.3) s we get 


h 


v 


u 


V 


C -/--/2 P - 


d VB 


(5-53 


The probability density function fg (h) of having a grain 

*=-<j | Ft j | j 

tip at any level h can be eb^ined~Tr©Pr--«quabioh~-4&r2d) -L ~ ' - 1 . 
by- substituting h in -place of z a -sinee -Thn&y- 


f R (h) 


(It + 1) (h + d ) k 

(2 k+1 - 1) df^ 5 

0 


0 <h < d 
- e 


(3 

(•5r6r) 


If M’ is the total number of grains lying below h , then, 
H f = M" Prob (h <h ) Gffe) 

U. 

where M n is the total number of grains lying below d^ 
for an average width of cut b. For a' nominal grain den- 


sity. 

C, M" 

is given by 



M n 

= -pr D 

b G 

% 

(-5V&-) 

Thus, 

M* 

= M n 

h u 

i frj (h) d h 

/ H 

05 • 93 


"0 
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Direct evaluation of h from the above equation 

is difficult, due to implicitness of the functions. There~ 

, ,1 ^ * 

fore, iterative method is used to solve equation 

in order to evaluate h for known values of k and d . 

"U. 0 

Once h is known for any cutting condition, can be 

\u, 

determined from equation (5v12d. 
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_ 5. %- - Surf ace Roughnes s. 

J r- l2 -i 

Farmer, Brecker/and Shaw j-7 1 ©*-; conducted tests 
under various grinding ''conditions and their experimental 
values of maximum peak to valley height are shown in 
Figure 5*3* It is observed that the maximum peak to 
valley height lie between the two dotted lines such that 
4 h <^h 6 h, where h is the centre line average value 

of the surface. Hence to a fairlj r good approximation, it 
can be assumed that 

h id ^ h. 


(5.1*0 
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Transverse distance , in 
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CHAPTER VI 


RESULTS, DISCUSSION AID CONCLUSIONS 
6.1 Introduction 

Dae system parameters k and d Q which control 
the probabilistic distribution of grains on the grinding 
wheel are the sources of variation for the random output 
of the grinding system (ground surface profile). Diere- 
fore, k and d have to be estimated so that the simu- 
lated and experimental profiles match. Die techniques 
adopted for evaluating these parameters are shown schema- 
tically in the block diagram (Figure 6.1). The values of 

d for various grinding conditions obtained from the 
6 

experimental profiles are fed to the variance analyser 
block. It is found (discussed later) that d g is practi- 
cally independent of the cutting conditions in the fine 
grinding range. With this value of d and a given 
value of k, the theoretical surface profile is simulated 
and fed to the comparator through correlator and spectral 
density analyser. The experimental surface profile is 
also fed to the comparator after being passed through 
correlator and spectral density analyser block. In the 
comparator block, the theoretical and experimental 
spectral estimates at different frequencies for various 
values of k are compared. Die outputs from the 
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comparator "block are then fed to the least square 
estimator where the values of k for any cutting condi- 
tion is estimated such that the comparator output is a 
minimum. This process is repeated for all cutting con- 
ditions. The least square estimates of k are then 
passed on to the variance analyser where the effect of 
d and v on k is evaluated. The regression estimates 
of k for its functional dependence on d and v are 
obtained in the regression analyser, the input for which 
are the least square estimates. The values of k, thus 
obtained, are then fed to the grinding system and the 
theoretical profile is generated which simulates the 
actual profile. The GLA value of the simulated surface 
profiles, can thus be calculated and compared with the 
experimental results. 

6.2 Experimental Details 

To obtain the experimental profiles, experiments 
have been carried out on a Hindustan Machine Tools Ltd. , 
Horizontal Spindle Surface Grinding Machine. The speci- 
fications of this machine are given in Appendix C. The 
following grinding conditions have been used: 

1 , Hhe el s A46-T5-V10. 

2. Workpiece: Mild Steel, length : 4 cm; width: 2 cm; 

height: 1 cm; hardness: Rg 60. 
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3* Wheel diameter s 30 cm 

b. Wheel depth of cut; 4- to 12 (microns) 

7. Table speed; 5 to 15 m/min. 

6. Wheel speed; 1500 r.p.m. 

7. Grinding fluid; dry 

8. Wheel dressing; (i) 2 passes at 12 microns and 

1 m/min cross-feed. 

(ii) 2 passes at 6 microns and 
1 m/min cross-feed 

(iii) spark out 

_2 

9. Nominal grain density; 160 cm 

10. Simulated width s 0.25 cm. 

All tests have been carried out under plunge cut 
conditions. The workpiece surface was mirror-finished 
and before performing each test ten passes were given 
with the dressed wheel on a seperate workpiece in order 
to eliminate the effect of "dulling" produced during 
dressing. This process was repeated for each test. The 
workpiece surface profile produced on the mirror- finished 
testpiece after a single pass was recorded on a Talysurf-4- 
Surface Measuring Instrument. As discussed in Chapter III, 
Sec. 3.2.1, a suitable 2x2 factorial experiment was 
adopted. A few representative traces of the experimental 
profiles are shown in Figures 3*3- • 



98 


6.3 Simulated Results 
6.3.1 Evaluation of d g 

To evaluate d and its possible dependence on 
d and v, the peak to valley height of the traces of 
the factorial experimental profiles were obtained. The 
values, thus obtained, are given in Table 6.1. 


TABLE 6.1 

Data for 2x2 factorial experiment for evaluating 

d e C/U ) 


Level of factor, d 

"i — r 

12 xtA. 

Level of factor, v 

5 5 m/min Si 5 m/rnin S5" 

m/min Si 5 m/min 

1 

Repetition 

2 

2.8 3.0 

3.0 2.9 

2.8 3.0 

3.0 3.2 

Total 

5.8 5.9 

5.8 6.2 

Grand Total 

23.7 



With the help of square identities and associated degrees 
of freedom, discussed in Chapter III, Sec. 3.2.4-, the ana- 
lysis of variance table is prepared (Table 6.2). 
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TABLE 6.2 



Analysis of 

variance 

Table for d^ 
e 


Source of 
Variation 

£> Sum of 

5 Sauares 

$ Degree of 5 Mean 

5 Freedom $ Squares 

5 Computed 
5 F-ratio 

Treatment 

combination 

(Tr) 

0.09+ 

3 

- 

- 

Factor, d 

0.0115 

1 

0.0115 

O .707 

Factor, v 

0.0315 

1 

0.0315 

1.877 

Interac tion 
(I) 

0.011 

1 

0.011 

0.677 

Within (S) 

0.065 

If 

0.01625 

- 

Total (T) 

0.119 

7 

- 

- 

From this table 

, it is observed that 

the computed F-ratio 


for d and v and their interactions are less than 7 . 71 , 
the value obtained from statistical table for ’F* distri- 
bution for upper 5 % significance level. Therefore it can 
be concluded that the effect of d and v on d 0 is negligibly 
small. Thus, the value of d can be taken to be equal to 
the mean of all values obtained from experimental profiles. 
Therefore, 

d e - £ ¥- 3 ' 0 / J - 
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6.3.2 Evaluation of k 

Once the value of d g is estimated, the following 
numerical values are calculated as 

N = Number of columns within the simulated width 
dr: 124, 

M = Number of grains per column 34, 
and the sample size, N and maximum lag number, are 

t alien as 500 and 50 , respectively. 

As shown in the block diagram representing the 
scheme of the solution technique (Figure 6.1), the theo- 
retical transverse surface profile is generated using 
known values of k and the random grain positions 
(equations (2.34) to (2.36)) following the computer flow 
charts (Figures- 4; ii and 2) . The output of the comparator block 
in terms of SSk (Chapter III, sec. 3*2>5)for iterative 
values of k is compared in the least square analyser 
for any grinding condition. The least square estimate of 
k is that value for which SSk is a minimum. This process 
is repeated for all cutting conditions. The comparison 
of spectral densities of theoretical profiles generated 
with the least square estimated values of k and expe- 
rimental profiles for some representative cases are shown 
in Figures 6.2. The least square estimates of k for 
the 2x2 factorial experiments are presented in Table 

6.3. 
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TAELE 6.3 

Data for 2 x 2 factorial experiment for evaluating k 


Level of factor,, d 5 

L 


1 

12 /u. 

Level of factor, ? i ^ 

m/min \ 15 m/mi nO 5 m/m±n$\5 m/mi n 


1 

5.0 

9.5 

8.0 

12.5 

Repetition 

2 

5.5 

8.5 

8.8 

13.0 

Total 

10.5 18.0 

16.8 

25.5 

Grand Total 



70.8 


The results of 

the analysis of variance are shown in 

Table 6 A. The procedure adopted 

was 

the same 

as that 

for evaluating 

d e* 

TABLE 6 A 




Analysis of 

variance table 

for k 


Source of 

5 Sum of 

IS Degrees 

of 5 

Mean 

<5 Computed 

Variation 

5 Squares 

5 Freedom 

5 

Sauares 

5 F- ratio 

Treatment 

combination 

5649' 

3 




(Tr) 





Factor, d 

23.805 

1 


23.805 

88.98 

Factor, v 

32.805 

1 


32.805 

122.70 

Interaction 





0,67 

(I) 

0.189 

1 


O.I 89 

Within (S) 

i .07 

If 


0.2675 

- 

Total (T) 

5746: 

7 


- 

- 
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This table shows that the computed F-ratio for 
d and v are much larger than 7.7 1, the value obtained 
from statistical table for 'F ! distribution for upper 
5 % significance level, where as for interaction effect' on 
k, the computed F-ratio is very small compared to the 
statistical value. Therefore, it can be concluded that 
the variation of k is due to the effect of d and v but not 
due to their interaction. . The functional dependence 
of k and d and v can be estimated by regression analysis. 

One of the requirements of the regression analy- 
sis, is the homogeneity of the variances in each trial. 

As discussed in Chapter III, Sec. 3.2.3, the variance of 
reproducibility of parameter, k is obtained and presen- 
ted in Table 6.7. 

TABLE 6.5 

Calculation of variance of reproducibility 


Trial 0 
Ilo . 5 

k i 

$ Repeated 

5 Trial No. 

5 k 5 

5 2 

5 k 

hi. 

1 

5.0 

5 

5.5 

5.25 

0.125 

3 

9.5 

7 

8.5 

9.0 

0.50 

2 

8.0 

6 

8.8 

8. if 

0.32 

if 

12.5 

8 

13.0 

12.75 

0.125 

S 2 (k) 

rep 

= 0 

.535 






2 

Since there is no variance, S^, (k) that 

2 

considerably exceeds the remaining ones, and (k), 

and S?, (k) are of the same order, homogeneity of vari- 
ances, can be assumed, however, to apply Bartlett test 
j~65 f, the value of c 2 ^- s calculated from equation (3*8) 
as 

Cg = 0.6182 

substituting the values of C 2 in Bartlett f s test sta- 
tistic, B, given by equation (3*7)? ue get 

B = 0.131 / 9<. 2 ( = 7.81) for 5% significance 

level and three degrees of freedom. Hence we can con- 
clude that the variances are homogeneous assuring the 
reproducibility of the values of k. 

Assuming multiple linear regression analysis, 
the equations (3.39), (3.40) and (3.41) given in Chapter 
III, Sec. 3.4, are solved for a-j, b^ and b 2 * The esti- 
mated regression of k on d and v is given by 

1 c = 1.35 + 0.412? d + 0.42 v (5.D 

^ere, a,, < and b 2 are 1.35, 0.4125 and 0.42, respec- 
tively. Solving equation (3.59), 997° confidence inter- 
val for and 3^ are 0.3087 <.A, <0.5163 and 0.35 < 
jS 2 0.618, respectively. 



6.3*3 Determination of Surface Roughness and 

Using equation (6.1) for various grinding 
conditions, the theoretical profile is generated and 
their spectral density curves are shown in Figures 6.3. 

From these theoretical profiles, the centre line average 
value (CLA) of the surface roughness are calculated. 

The CLA values, thus obtained, for various cutting con- 
ditions for "both theoretical and experimental profiles 
are shown in Figures 6.4- and 6.5* Using the computer 
flow charts,, the active number of grain density respon- 
sible for surface generation are numerically obtained and 
plotted in Figures 6.6 for various grinding conditions. 

6.4- Results of Approximate Analysis 

Equation (5-13) of Chapter V, Sec. 5.2, is solved 
iteratively with an initial value of h^ as 0.5/^- with 
convergence criterion 0.001 . The values of centre- 

line average (CLA) of the surface roughness are calcula- 
ted by using equation (5»l4-). The number of active grain 
density, C^, can be obtained from equation (5.12). The 
variation of surface roughness with d and v are shown 
in Figures 6.7 while the variation of active grain density 
with h u is shown in Figures 6.8 
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6.5 Discussion 

As the surface generated is the sample realisation 
of grains on the wheel surface, the surface profile produ- 
ced during the wheel work interaction will he random in 
nature. It has been shown through statistical analysis 
that d is independent of the cutting condition, there- 

G 

fore, this randomness is primarily due to the effect of 
radial distribution parameter. The number of cutting 
points on the wheel surface though- varies with the radial 
distance from the outermost grain, it quickly gets satu- 
rated as - the radial depth is increased. Therefore, it is 
quite reasonable that the effective profile depth remains 
more or less constant for a particular wheel and a large 
part of fluctuations in the experimental data may be due 
to chance variation. 

When the wheel is stationary, the depth of cut 
and the table speed are zero and the value of k is a^ 

• which shows the static distribution of the grains on 
wheel surface . Earlier analysis seve the value of 

k under the above circumstance as zero which does not 
have any physical bearing. The dynamic influence of table 
speed and depth of cut on radial distribution parameter 
is functionally expressed by equation (6.1). 

Actual ground surface is distorted from an ideal 
profile by various causes, some of which are probabilistic 
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and others are systematic. Comparison of real surface 
with ideal one will provide information about the real 
system if the ideal one is properly simulated to satisfy 
the spectrum characteristics associated with the real 
surface. In earlier analyse s , j 57 ? 61 j, it was observed 
that the power spectra of the actually ground surface is 
very much amplified in low frequency range having peaks 
at very small frequency -and with increasing values of d 
and v, the profile variances increased. In the results 
shown in Figures 6.2 the above criterion is observed for 
the experimental profiles, having nigh spectral density 
in low frequency range (less than 0.1 5 c/Zi ) with peaks 
concentrated around 0.02 c/ A . The power spectral den- 
sity of the simulated profile shows amplification in the 
low frequency range (less than 0.1 5 o//\ ) but no peaks 
are observed (Figures 6.3) • 

A possible explanation for this, lies in the fact 
that the sample experimental surface, even though mirror 
finished, to start with, will have some irregularities 
and is not ideally flat as considered for the theoretical 
analysis. The effect of vibration during grinding and 
individual grain deflection will have contribution to the 
above effect. Therefore, the error involved in estimation 
of the power spectral density of the real system in terms 
of the theoretical spectral density will be 
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C(fq) = 'G (fq, ^ ' G (fq, £ , k) (6.2) 

Our aim is to obtain a sound, theoretical model, so that 
this error is vexy small. However, one might want , for 
example, to minimise the probability 

Prob < I'S- (fq, ^ ) - Q' (fq, f , k) I at a<<1 

l ' - J \ y y 

( 6 . 3 ) 

Such a requirement, though very realistic, is difficult 
to handle analytically. In the model, this criterion is 
reasonabljr solved by the minimisation of the mean square 
error, e, j 62 given by 

e = E jj'G (fq, ^ ) - G (fq, > # k) j ~] (6 A) 

■where it was assumed that £(fq) is uncorrelated. Least 
square estimate of k is obtained from the above equation. 

As the depth of cut, d and v are increased, the 
profile variances represented by the area under the curve 
increases (^Figures 6. ; 3 ). The power spectral density also 
increases. As the root mean square value of the profile 
is equal to the area under the spectral density curve, 
the increase in profile variance indicates relatively rough 
surface with increasing d and v. 

The effect of grinding conditions on surface 
roughness is shown in Figures 6Jf and 6.*?. It is seen 
that with the increase in the value of d and v, surface 
finish deteriorates. The two important parameters that 
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will have significant effect on surface roughness are the 
undeformed chip thickness and active grain density. With 
increasing values of chip thickness, the peak to valley 
height of surface profile will increase while the increase 
in grain density will tend to decrease it. It appears 
that the chip thickness has more pronounced effect causing 
a deterioration of surface finish. This aspect is indica- 
ted by both theoretical and experimental results shown 
in the above figures. Further, there appears to be a 
reasonably close agreement between the experimental and 
simulated results. 

The effect of grinding conditions on active grain 
density is shown in Figures 6.6a and 6.6b. These figures 
clearly indicate that both d and v affect the number of 
cutting points responsible for surface generation. An 
individual grain will only see the undeformed chip thick- 
ness and this value will increase with increasing values 
of d and v. As the chip thickness increases, more and 
more grains will get involved in cutting. This, however, 
will not increase infinitely but will tend to saturate 
beyond a certain value . This aspect Is clearly seen in 
the above figures and the trend of the result agrees with 
those published earlier 

Once the dependence of k on d and v are known, 
it is possible to develop an approximate analysis for 
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evaluating the uncle fo rued chip thickness, t. In deriving 
this expression for t, it has been assumed that there 
exists an average linear grain spacing, /\ . However, the 
radial positions of grain have been assumed to be distri- 
buted randomly between 0 and t. The upper limit for t 
has been taken to be that value where G gets saturated. 
Thus, the effective profile depth has been assumed to be 
the limiting value of t for the approximate analysis. 

In previous analyses {33., 28, 29 j , an average linear 
grain spacing has been taken while all grains have been 
assumed to be at the same radial height from the wheel 
surface. In addition, it was also assumed that the actual 
grain density C remains constant and independent of t. 


Results of the approximate analyst sCFigures 6.7) 
show that the surface finish deteriorates as the depth of 
cut and table speed are increased. Comparison with expe- 
rimental results also show reasonable agreement* The 
plot of C versus t ( Figure 6.8) clearly indicates that 
the active grain density increases with increasing value 
of t and agrees with the trend of the experimental curve 
obtained by Nakayama and Shaw j 33 J* Such a curve is 
valid only in the fine grinding range since the active 
grain density can not rise infinitely but must saturate 
at some value of t. Therefore, the approximate analysis 
suggested is valid only in the fine grinding range. 
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6 .6 Conclusions 

A probabilistic approach seems to be natural 
for forming an analytical basis for the grinding process 
which involves several complex random parameters. 

The analysis shows that the process can be des- 
cribed in terms of two fundamental parameters, the radial 
distribution parameter and the effective profile depth. 
Results show that the radial distribution parameter k 
is a linear function of the wheel depth of cut and table 
speed, while the effective profile depth is statistically 
independent of cutting conditions. The computer simulated 
surface profiles show close agreement with the experimental 
profiles in terms of stochastic characteristic such as 
spectral density functions. Results also indicate that 
surface finish deteriorates with increasing values of 
wheel depth of cut and table speed and agrees with the 
experimental results.. 

Once the dependence of radial distribution para- 
meter on cutting conditions is established, the approximate 
analysis proposed can be used to predict the values of 
chip thickness, and surface finish. Comparison of theo- 
retical and experimental results show reasonable agreement. 

In the simulated analysis, only single pass 
plunge cut conditions have been considered and also the 
workpiece has been assumed to be a perfectly smooth 
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surface. This analysis can he extended to include the 
effect of original workpiece roughness, on the surface 
generated. In single pass operations the wear of abra- 
sive grains can he neglected hut in multipass operations 
this effect must he incorporated. In all fine grinding 
operations, spark out condition are used during the last 
few passes. The work can be further extended to include 
these effects. 

The mechanics of chip formation with single 
grains is still not clearly understood. It is believed 
that there are three regimes, rubbing, ploughing, and 
cutting, during cutting with single abrasive grains 
depending upon the depth of cut. In the present analysis, 
it has been assumed that every interference between the 
grain and the workpiece removes material in the form of 
a chip without side pile-up. Once this mechanism is 
clearly understood, this aspect can also be incorporated 
in the analysis. 
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APPENDIX A 


The random variables p^, p^, . etc. are 
constructed with the density function given by equation 
(2.$). Assuming that the random variables are indepen- 
dent, the distribution function of p is obtained from 

P_ 

F (p) = [' — § (L - p) d p (A.1) 

P V L 

o 

where L is the maximum width ox cut. If, Rv. , is a 

3 

random number, uniformly distributed between 0 and 1 , we 
have ! 69 | 



F p (P) 



(A. 2) 


Solving the equation (A. 2), we get 


p = L (1 ±\J 1 - R u _) (A. 3> 

3 

As the maximum value of axial pitch will be less than or 
equal to L, 

p = L Cl -s/1 - R " ) (A.4) 

0 

Hence, by definition of axial pitch, 

x. = x • j , + L ( 1 - VI - R ) (A. 5) 

3 3“ * 
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The distribution function, F y ( /j.y), can be 

obtained from equation (2*31). Starting from an arbi- 
i.«~i 

trary point y = 0, we form the points 





(A. 6) 



Ay 1 = - ?l 1 log e (1 - \) 

= - A. log (R. ) <A.9) 

i e u jL-1 

where, 

R = 1 - R (A. 10) 

" 1-1 U 

is also random number uniformly distributed between 0 and 
1. Hence, the random point given by equation (A.6) can 
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be obtained using equation (A. 10) as 


y 1 = y 1_1 -A 1 log (R ) 

' e i-1 


(A. 11 ) 


All other points can be similarly obtained. Using the 

. , i i+1 M 

same principle , the random points , a , z , . . . , z can 

be generated from the equation 

i 




z 

= i f z (z 1 ) d z 
o 


(A. 12) 


which on simplification gives 


l/(k+1 ) 


A = A [ (2 k+1 - 1) B + l") 


e V 


T i t 

u i J 


e 

(A.13) 


for i = 1, . M. 



APPENDIX E 


In Figure E.1, E (i, i) is the intersection of 
the path traced by ith and jth grain in any column. The 
coordinates of the lowest point of ith grain in its path 
is given by (y 1 , z 1 ) in the plane along the cutting direc- 
tion. Similarly, the coordinates of the lowest point of 
jth grain in its path will be (y%z^). As shown in Figure 
B.1, A1 and A2 represent the radii of the locus of paths 
of ith and jth grain respectively. Therefore 

A1 2 - yl 2 = A2 2 - y2 2 (B.1) 

or 

y2 2 - yl 2 = A2 2 - Al 2 (B.2) 

viie re, 

Al = R - z 1 (B.3) 

and 

A2 = R - z\ (B.1+) 

and 

yl + y2 = y 3 - y 1 (B.5) 

Solving equations (B.2) and (B.5), we get 

C7 3 - r 1 ) - ^ (B - 6) 

- r_J 

Again, 

A2 2 - Al 2 = (R - z^ ) 2 - (R - z 1 ) 2 

■ ■ . 2 

= - 2R (z*^ - z 1 ) * (z' 3 " - z 1 ) 

(B.7) 




-2 


Neglecting the higher order terra (z J - z 1 ) in equation 
(B.7), we get 

A2 2 - Al 2 ^cr- 2R (z^ - z X ) ( 2 . 8 ) 

Therefore, 


yi 


1 [ > C . z.. J .= |7 - )- + (y 3 _ y i) j (B.9) 


2 L (y^ - y 1 ) — 

and the y„ coordinate of the intersection point in the 
(y, z) plane will "be 


, . - \ i , 1 2R (z^ - z X ) / j _J 

y (i, 3 ) = y + 2 — 71 rr Cy “ y 

^ (y d - y ) 


7 


(B. 10) 

The vertical height of the intersection point will be 


z (i, j) = z + A3 

\diere A3 can be obtained from equation 


(B.11) 


i\2 


(R - z x - A3) 2 + yi 2 = (£!.- z 1 ) 

which on simplification reduces to 

,2 


(B.12) 


A3 t 

Therefore, 

z (i, 3 ) 


_Xt 


2 (R - z 1 ) 


(B.13) 


l . 
Z 4* 


(y (i»3) - y 1 ) 2 


(B.14-) 


2 (R - 

Using the same procedure, the coordinates of the inter- 
section of grooves in transverse section at any distance 

y can be obtained, where the circular shaped grooves of 
■ o * 
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radius f' are superimposed with its mid point at the 
lowest point of the groove in any column. Therefore, 
it can be shown that 


*1,3 “ X i + 2 


1 I g (a i - Z 1 ty n » 


Cij - xTT 


+ (x i - x fJ 

(B.15) 


and 


1, J 


= z- (y) + 


(x i.i - x i } 


i •'o' 


g 


(B.16) 





APPENDIX 


n. 




MACHINE SPECIFICATIONS 




Range Value 

Units 

1. 

Working Surface 

250 x 1000 

mm 

2. 

Maximum Transverse Traverse 




of Grinding Wheel 

320 

mm 

3. 

Maximum Longitudinal Traverse 




of Table 

1150 

mm 

4. 

Maximum Grinding Length 

1000 

mm 

5. 

Maximum Grinding Height 

400 

mm 

6. 

Infinitely Variable Table Speeds 

2.5 - 25 

m/mi n 

7. 

Automatic Continuous Traverse 




Feed of Grinding Wheel 

0.3 - 3.5 

m/min 

8. 

Automatic Intermittent 




Tranrver.se Feed of Grinding Wheel/ 




Table Stroke 

1 - 32 

mm 

9. 

Downfeed of Grinding Wheel 




1 Graduation = 

0.002 

mm 

10. 

Rapid Vertical Speed of Grinding 




Wheel 

0.6 

m/min 

ii. 

Grinding Wheel Speed 

1440/2880 

rpm 

12. 

Wheel Diameter Upto 

300 

mm 

13. 

Wheel Width Upto 

63 

mm 

lb. 

Wheel Bore 

76 

mm 

15. 

H.P. of Grinding Motor 

4/5.7 

h*p* 



